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Abstract 


AT20211fa, also known as ZTF21aayokph, was detected by the Zwicky Transient Facility on 2021 May 4, at 
05:34:48 UTC. Follow-up observations were conducted using a range of ground-based optical telescopes, as well 
as Swift/XRT and VLA instruments. AT20211fa is classified as an “orphan afterglow” candidate due to its rapid 
flux decline and its reddened color (g — r = 0.17 + 0.14 mag). For an optical transient source without prompt 
gamma-ray detection, one key point is to determine its burst time. Here we measure the burst time through fitting 
the initial bump feature of AT2021lfa and obtain its burst time as 2021 May 3, at 22:09:50 UTC. Using 
afterglowpy, we model the multi-band afterglow of AT2021l1fa and find that the standard model cannot 
reproduce the late radio observations well. Considering that the microphysical parameters €e, €g (the energy 
fraction given to electrons and magnetic field), and £y (the fraction of accelerated electrons) may vary with time, 
we then model the afterglow of AT20211fa taking into account the temporal evolution of the physical parameters 
€e, €g, and £y and find in this case the multi-wavelength observations can be reproduced well. The initial Lorentz 
factor of AT20211fa can be estimated from the peak time of the early afterglow, which yields a value of about 18, 
suggesting that AT20211fa should be classified as a “dirty fireball.” From the upper limit for the prompt emission 
energy of AT20211fa, we obtain that the radiation efficiency is less than 0.02%, which is much smaller than that of 
ordinary gamma-ray bursts (GRBs). It is also interesting that the fitted values of jet angle and viewing angle are 
very large, 0. ~ 0.66 rad, 0, ~ 0.53 rad, which may lead to the low Lorentz factor and radiation efficiency. When 
compared with GRB afterglow samples, it is evident that the onset bump timescale of AT20211fa satisfies the 
empirical relationships observed in GRB samples. Additionally, the luminosity of AT20211fa falls within the range 
of observations for GRB samples; however, approximately 1 day after the burst, its luminosity exceeds that of the 


majority of GRB samples. 
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1. Introduction 


Gamma-ray bursts (GRBs) represent the most energetic and 
brightest explosive events in the universe, occurring at cosmo- 
logical distances. GRBs can be classified into two categories 
based on the duration of gamma-ray emission: long gamma-ray 
bursts (LGRBs), which have durations greater than 2 s, and short 
gamma-ray bursts (SGRBs), which have durations less than 2 s 
(Kouveliotou et al. 1993). Typically, LGRBs originate from the 
core-collapse of massive stars, whereas SGRBs arise from the 
merger of compact binary systems, such as neutron star-neutron 
star systems or neutron star-black hole systems (Eichler et al. 
1989; Woosley 1993; Zhang & Mészáros 2004; Kumar & 
Zhang 2015). In the classical GRB fireball model, a central 
engine, formed by the collapse of a massive star or the merger of 
compact objects, launches a collimated and ultra-relativistic 
outflow (Paczynski 1986; Piran 1999; Racusin et al. 2009; Zhang 
et al. 2015). Within this outflow, dissipative processes result in an 
intense prompt gamma-ray emission lasting from seconds to 


minutes (Rees & Meszaros 1992; Thompson 1994; Maxham & 
Zhang 2009; Hu et al. 2014). Furthermore, interactions between 
the outflow and the surrounding circumburst medium produce a 
long-lasting, broadband synchrotron radiation known as afterglow 
emission, persisting from days to months (Mészáros & Rees 1997; 
Sari & Mészáros 2000). 

With the launch of GRB satellites equipped with rapid 
response capabilities, such as Swift (Gehrels et al. 2004) and 
Fermi (Meegan et al. 2009), an increasing number of GRB 
afterglows have been detected. These afterglows are identified 
through follow-up observations triggered by gamma-ray 
emissions. However, there are instances of "orphan after- 
glows," where GRB afterglow radiation is detected without a 
corresponding prompt gamma-ray emission. This phenomenon 
can be attributed to two primary factors. First, GRB jets exhibit 
a high degree of collimation, making it challenging for GRB 
satellites to observe the prompt emission when the viewing 
angle exceeds the half-width angle of the jet core, a scenario 
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referred to as the “off-axis” case. However, during the 
afterglow phase, as the bulk Lorentz factor decreases, the 
radiation's beaming angle widens, which enables the observa- 
tion of afterglow radiation to become possible (Huang et al. 
2002). Second, “dirty fireballs,” characterized by initial Lorentz 
factors much less than 100 due to significant baryon loading, 
can produce broadband afterglow radiation. However, the 
associated prompt emission is too weak for detection by GRB 
satellites (Paczyński 1998; Dermer et al. 1999). Consequently, 
in order to understand the properties of the whole GRB 
population well, it is essential to search and study the properties 
of the “orphan afterglows." 

Before the deployment of the Zwicky Transient Facility 
(ZTF; Bellm et al. 2018; Graham et al. 2019), only three optical 
afterglows had been observed which were not triggered by 
GRB prompt emission. In a retrospective survey, the Odyssey, 
INTEGRAL (SPI-ACS), and Konus missions identified GRB 
140226A as the prompt emission counterpart to iPTF1l4yb 
(Hurley et al. 2014; Cenko et al. 2015), while POLAR, 
AstroSat, Konus-Wind, and INTEGRAL observed that GRB 
170105A was the prompt emission counterpart of ATLAS17- 
aeu (Bhalerao et al. 2017; Svinkin et al. 2017). PTF11agg, the 
only optical transient source without a GRB counterpart (Wang 
& Dai 2013; Wu et al. 2014; Wang et al. 2015, 2016), was 
identified as likely being a “dirty fireball” by Cenko et al. 
(2013). The ZTF, characterized by its high-cadence, large field 
of view, and rapid follow-up capabilities, has enabled people to 
find more optical afterglow candidates (Dekany et al. 2020). 
From 2020 to 2022, ZTF detected twelve optical afterglow 
candidates, four of which (AT2020blt, AT2021any, AT21Ifa, 
and AT2022cmc) did not exhibit corresponding GRB prompt 
emissions in retrospective searches. Regarding AT2020blt, 
Sarin et al. (2022) proposed that it might be interpreted as the 
on-axis afterglow of a long-duration GRB. Due to its low 
gamma-ray radiation efficiency, falling below 0.396—4.596, its 
prompt emission cannot be detected by current missions. The 
multi-band afterglow data fitting conducted by Xu et al. (2023) 
for AT2021any considered both the top-hat jet model and the 
structured Gaussian jet model. Their analysis revealed that the 
structured Gaussian jet model provided a better fit result, 
suggesting that AT2021any might be a “dirty fireball” observed 
on-axis. Furthermore, Andreoni et al. (2022) confirmed that 
AT2022emc is a jetted tidal disruption event based on the 
characteristics of its X-ray, submillimeter, and radio light 
curves. 

The light curve of AT2021lfa exhibits features consistent 
with those of GRB afterglows. However, the absence of 
corresponding prompt emissions suggests that it may be a 
potential candidate for an “orphan afterglow.” Lipunov et al. 
(2022) employed the smooth optical self-similar emission 
(SOSS) model for GRBs to model the optical light curve of 
A'T20211fa. The fitting results revealed an initial Lorentz factor 
of T'o = 20 + 10. Consequently, they proposed that AT20211fa 
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is most likely to be a “dirty fireball.” Ho et al. (2022) utilized a 
single power-law function to fit the light curves of the r, g, and 
i bands of AT20211fa. 

In this paper we re-investigate the properties of AT20211fa 
by incorporating more observational data from radio, optical 
and X-ray bands, thereby imposing further constraints on the 
nature of AT20211fa. Based on the empirical formula for the 
early-phase afterglow bump, we determine both the burst time 
and peak time of AT20211fa. Employing the afterglowpy 
tool, we conduct a multi-band afterglow fitting analysis for 
AT20211fa. Throughout the fitting process, we consider the 
temporal evolution of microphysical parameters e, €g (the 
energy fraction given to electrons and magnetic field), and Ey 
(the fraction of accelerated electrons) which significantly 
improves our fitting results, particularly in the radio band. 
The fitting results strongly suggest that AT20211fa is highly 
likely to be a “dirty fireball” with an initial Lorentz factor of 
approximately 18. In Section 2 we describe the multi- 
wavelength observation of AT2021lfa. The determination of 
A'T20211fa's burst time and the multi-band light curve and 
spectra fitting are presented in Section 3. In Section 4, we 
discuss the property of AT20211fa and compare its light curve 
features with those of GRB afterglows. Finally, we provide a 
summary of our findings and conclusions in Section 5. 


2. Observations 


AT20211fa (ZTF21aayokph) was discovered by ZTF at the 
position R.A. (J2000): 12:32:48.725, decl. (J2000): —01:29:22.56. 
The optical transient was initially detected at 05:34:48 UTC on 
2021 May 4, with a magnitude of r= 18.60 + 0.08, and at 
07:30:00 UTC on 2021 May 4, with a magnitude of 
g= 18.80 +0.11. The most recent upper limits were recorded 
on 2021 May 2, at 07:27:50 UTC (referred to as fup), with 
r > 20.23 mag, and at 06:09:56 UTC with g > 20.58 mag (Yao 
et al. 2021b). The Liverpool Telescope (LT) (Steele et al. 2004) 
conducted rapid follow-up observations in the r filter and g filter at 
15.5 and 19 hr after the first detection, revealing a rapid fading of 
AT2021 Ifa (1.9 mag day! in r-band) and colors (g — r= 0.17 + 
0.14 mag) consistent with a non-thermal spectral energy distribu- 
tion (SED). Considering the fast decay of flux and the reddened 
color, Yao et al. (2021b) suggested that AT20211fa is a potential 
candidate for a GRB optical afterglow. 

In the optical band, in addition to the observations by ZTF 
and LT, DDOTI (Watson et al. 2021), RATIR (Butler et al. 
2021), NEXT (Fu et al. 2021), Assy (Kim et al. 2021), Lowell 
Discovery Telescope (O'Connor et al. 2021), Mondy and 
Koshka (Pankov et al. 2021), and LBT (Rossi et al. 2021) also 
conducted rapid follow-up photometric observations of 
A'T20211fa. Of particular significance, the MASTER-OAFA 
robotic telescope had observed AT2021lfa in survey mode 
three hours before the first detection by ZTF (Lipunov et al. 
2021). MASTER captured the phase of steady brightening in 
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AT20211fa, providing valuable insights for constraining the 
time of peak luminosity. Yao et al. (2021a) observed 
A'T20211fa with GMOS-S and identified z = 1.063 as the most 
probable redshift. Optical photometry has been corrected for 
Galactic extinction with E(B — V) — 0.026 mag (Schlafly & 
Finkbeiner 2011). The amount of extinction Ay of the host 
galaxy, obtained by fitting the SED of AT2021lfa with the 
SMC template extinction law (Li et al. 2008), is very small and 
can be neglected. 

In the X-ray band, Swift/XRT observed AT20211fa 1.0 day 
after the initial ZTF detection, with two separate epochs 
totaling an exposure time of 5ks. X-ray emission from 
AT20211fa was detected only in the first epoch. For observation 
data, we refer to Ho et al. (2022). Assuming the photon index 
To = 2, we converted the observed flux in the 0.3-10 keV range 
to the flux density at 2 keV. 

In the radio band, VLA observed AT20211fa at nine different 
epochs. The observations spanned from a few days after the 
first ZTF observation to several hundred days later. We 
obtained the radio band data from Ho et al. (2022). 


3. Method and Results 
3.1. Fitting the Burst Time of AT2021lfa 


Precisely determining the burst times (Tọ) is imperative for 
investigating the fundamental physical mechanisms responsible 
for both the prompt gamma-ray emission and the afterglow 
radiation of GRBs. Traditionally, the burst time is determined 
by the trigger time observed from high energy satellites. 
However, in cases involving “orphan afterglow,” confirming 
the burst times presents a challenge. Ho et al. (2022) utilized a 
single power-law function to fit the light curves of the r, g, and 
i bands of AT20211fa. This allowed them to estimate the burst 
time as 79, — 2021 May 3, 06:57:36 (all dates and times 
provided in this article are in Coordinated Universal Time, 
UTC). In contrast, Lipunov et al. (2022) employed the SOSS 
model to match the light curve of AT2021lfa and derived a 
burst time: Twig = 2021 May 4, 01:33. A discrepancy of 
0.77 day is evident in the burst times obtained through these 
two approaches, which would significantly impact the 
subsequent analysis of light curve fitting for AT20211fa. 

To further determine the burst time of AT2021lfa, we 
consider an alternative method for estimating this parameter. In 
the framework of the GRB internal and external shock model, 
when the fireball reaches the deceleration radius, it undergoes a 
transition from the coasting phase to the deceleration phase. 
During the early afterglow phase, the jet interacts with the 
surrounding circumburst medium, leading to electron accelera- 
tion and then producing synchrotron radiation. Consequently, 
the afterglow light curve begins to rise. Once approximately 
half of the fireball’s energy is transferred to the circumburst 
medium, the light curve reaches its peak and subsequently 
enters a declining phase. As a result, the light curve would 
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Figure 1. The light curve of AT20211fa in r band. The solid line is the fitted 
result with the empirical formula proposed by Kocevski & Liang (2001). 


exhibit a gradual bump during the initial afterglow phase. For 
A'T20211fa, a noticeable bump is observed. Kocevski & Liang 
(2001) introduced an empirical formula to model the early light 
curve bump of the afterglow 


r r+1 = Fi 
F()- E, t — dt d ! r t — dt 
tyi — dt r+d r+ d\tp — dt 


(1) 


where F, represents the maximum observed flux during the 
peak time, and ¢,, is defined as the time span between tup and 
the peak time fpeax, calculated as t5 = tpeak — fup. The temporal 
gap from fup to the burst time Tọ is denoted as dt, dt = To — fup. 
The parameter r represents the rising power-law index, while d 
signifies the decaying power-law index. 

The magnitude of AT20211fa, observed by the MASTER 
telescope, is calibrated by using reference stars within the G 
band from the Gaia Early Data Release 3 (EDR3) a 
wavelength range similar to that of the r band. Therefore, we 
approximate the observations in the G band as those in the r 
band. Furthermore, we convert the initial observational data 
from the ZTF g band to the r band by taking the spectral index 
as — 1 (the spectral index of AT20211fa was measured to be —1 
at about 04:00 UTC on May 5). Subsequently, we employ 
emcee, a Python-based ensemble sampling technique for 
Markov Chain Monte Carlo, to fit the r-band light curve of 
AT20211fa using Equation (1) by taking F,, tjj, dt, r, and d as 
free parameters. 

The fitted light curve is displayed in Figure 1, while Figure 2 
illustrates the posterior distribution of parameters. From 
Figure 1, it can be observed that after about 4 x 10? s, the light 
curve significantly flattens, possibly due to contributions from 
the host galaxy. Therefore, observational data beyond 4 x 10^ s 
were not included in our fitting. The obtained parameters are as 
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Figure 2. The posterior parameter distribution for AT20211fa. 


follows: the peak flux F, = (1.560*0055) x 107^ Jy; the dura- 
tion £j; between the reference time fup and the peak time fpeak is 
(1717317655) s; the time interval dt between the burst time Tp 
and the reference time fup is (13931672235) s; the rising power- 
law index r = (3.007*1139); and the decaying power-law index 
d = (1.95470 6%"). Therefore, the burst time is determined as 
Ty = 2021 May 3, 22:09:50797:90:33, and the peak time is 
represented by fpeak = 2021 May 4, 07:10:01700:11:32. Compar- 
ing our results with those obtained by Lipunov et al. (2022), it is 
evident that our inferred burst time is earlier than theirs and 
meanwhile the peak time is somewhat later. The fitted rising 
power-law index, r~3, is consistent with the rising index 
expected for the early afterglow light curves in a uniform 
medium (Panaitescu & Vestrand 2008). 


3.2. Fitting of AT2021lfa Multi-band Afterglow 


Afterglowpy is an open-source modeling tool designed 
for GRB afterglows. It enables simultaneous fitting of multi- 
band data in X-ray, optical, and radio frequencies (Cunningham 
et al. 2020; Ryan et al. 2020; Li et al. 2022). The tool employs 
the single-shell approximation to simulate the interaction 


between GRB jets and a uniform external medium. Addition- 
ally, it utilizes a semi-analytical approach to calculate the 
synchrotron radiation generated by the forward shock. After- 
glowpy incorporates not only the top-hat model but also the 
Gaussian jet model and the Power-Law jet model. 

The fluxDensity function in afterglowpy incorpo- 
rates ten independent physical parameters essential for better 
fitting of the afterglow light curve. These parameters include 0, 
(viewing angle), Eo (on-axis isotropic equivalent energy), 6. 
(half-width of the jet core), 0,, (truncation angle of the jet’s 
"wing"), b (power for power-law structure), n (circumburst 
density), p (power-law index for electron distribution), e, 
(thermal energy fraction in electrons), €g (thermal energy 
fraction in the magnetic field), and £y (fraction of accelerated 
electrons). We employ the emcee Python package, in 
combination with the fluxDensity function from after- 
glowpy to fit the afterglow light curve. We assume log- 
uniform prior distributions for Eo, n, €e, and eg, with the 
remaining parameters assigned uniform prior distributions. 

The standard afterglow model assumes that microphysical 
parameters, such as €,, eg, and £y, remain constant throughout 
the entire afterglow evolution. However, studies have indicated 
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Figure 3. Multi-band afterglow data and fitting of AT20211fa, where squares represent observed data points not utilized for fitting, and downward triangles denote 


upper limits. 


that these parameters may vary with time during the course of 
the afterglow's evolution (Yost et al. 2003; Fan & Piran 2006; 
Ioka et al. 2006; Wei & Fan 2007). By fitting the early 
afterglows of several GRBs, it has been observed that the 
microphysical parameters e, and eg of the ultra-relativistic 
forward shock and mildly relativistic reverse shock are different 
(Fan et al. 2002; Kumar & Panaitescu 2003; Panaitescu & 
Kumar 2004; Wei et al. 2006), which suggests that these 
parameters may depend on the strength of the shock. To take 
into account the time-dependent variations of €e, €g, and £y, our 
fitting procedure incorporates the assumption that these three 
physical parameters change over time following power-law 
behaviors: e, = €,ot^,, eg = egot^?, and £y = Eyot™, where Eeo, 
€go, and £yo represent the initial values of the respective 
parameters, while a), a2, and a3 signify the corresponding 
power-law exponents. 

We employed afterglowpy to fit observational data in the 
radio, optical, and X-ray bands for AT2021Ifa after the peak 
time. We considered three jet models for data fitting: top-hat, 
Gaussian, and Power-Law, respectively. As illustrated in 
Figure 3, the observed flux in the r and i bands exhibits a 
flattening trend beyond 4 x 10? s which may be the contribution 
from the host galaxy, thus we excluded this segment from the 
fiting process. It is worth noting that the intricate fluctuations 
observed in the radio band prior to 12 days may be attributed to 
interstellar scintillation (Rickett et al. 1984), therefore these radio 
data were excluded from the fitting analysis. The results of the 


Table 1 
The Results of Fitted Physical Parameters for AT20211fa 
Model Top-hat Gaussian Power-law 
6, [rad] 0.526 0181 0.521*0.181 0.530*015 
log, Eo 54.714043 54.928*0496 54.848*0406 
6, [rad] 0.662 0207 0.905*0206 0.907020 
0, [rad] A: 0.69770 312 0.687 0208 
b x B 5.055333 
log,ono 1.039* 04627 1.139035 1.04170 360 
p 3.092003 3.091*0083 3.085*0055 
ai —0.055 0.059 —0.060 5075 —0.058* 0074 
logio co —1.17 749318 ~ 1.170793 = 1.115402) 
a —0.355* 0065 —0.390* 0055 —0.382* 0,070 
logio €B0 —4465*0 5 —4.434 0407 —4402 515 
az —0.204*008 —0.220*0024 —0.207*0030 
Eno 0.695 0353 0.7297 533) 0.728185 


fitting procedure are summarized in Table 1. Remarkably, in 
both the Gaussian and Power-Law jet models, the half-width of 
the jet core exceeded the truncation angle of the jet, which 
indicates that the jet structure should be uniform. The physical 
parameters derived from modeling using the top-hat jet model 
are as follows: viewing angle 0, = 0.52670 181 radians, isotropic 
equivalent energy log), Eo = 54.7741 6433 , half-width of the jet 


core 0. = 0.662702? radians, power-law index for electron 
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Figure 4. The posterior distribution of physical parameters of the top-hat jet model for AT20211fa. 


distribution p = 3.092'0024. circumburst density logio no = 
1.039+9:637 , the thermal energy fraction in electrons and the 


magnetic field are logi, e,9 = ed. 77020 and log, epo = 


—4.46510- 107 respectively, the fraction of electrons that get 


accelerated Eyo = 0.695*02/5, and their corresponding temporal 


evolution indexes are aj = —0.055^0059. a, = —0.355 005. 


az = —0.2047 003. respectively. 

The multi-band light curve of AT20211fa was successfully 
fitted using the physical parameters derived from the top-hat jet 
model, as presented in Figure 3. Notably, the wu-band 
observations exhibit a lower trend compared to the fitted light 


curve, possibly due to Lya absorption. Furthermore, Figure 4 


illustrates the posterior distributions of the fitted parameters. 
Figure 5 illustrates the SED of AT20211fa at 1.24 days. The 
SED was modeled using afterglowpy. It is evident from the 
figure that the SED exhibits a break at approximately 
4 x 10/5 Hz (the cooling frequency v,). Consequently, the 
optical and X-ray observations fall within the v, < Von < 
Vy < Ve region at 1.24 days. 


4. Discussion 


When a transient source exhibits similar light curve features 
to those of a GRB afterglow, but no associated prompt gamma- 
ray emission is detected, it is categorized as an “orphan 
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Figure 5. Observed SED of AT20211fa (points) and the top-hat fit (red line) 
from afterglowpy at 1.24 days. The dotted line indicates the cooling 
frequency (ve). 


afterglow.” The “orphan afterglow’ phenomenon can be 
attributed to two potential factors. First, these events could 
arise from what are known as “dirty fireballs.” Second, they 
might result from off-axis observations, where the GRB jet is 
not directed toward Earth. 

The Third Interplanetary Network (IPN3) comprises ten 
space satellites designed for continuous all-sky monitoring of 
GRBs. The most sensitive GRB satellites within IPN3 include 
the KONUS instrument on the Wind spacecraft, the Swift Burst 
Alert Telescope, and the Fermi Gamma-ray Burst Monitor. For 
AT20211fa, Ho et al. (2022) performed a search within the 
IPN3 GRB observation archives to identify potential prompt 
emission counterparts. However, no corresponding prompt 
emission was detected. As a result, AT2021lfa has been 
classified as an "orphan afterglow" candidate. Based on 
KONUS-Wind observations, Ho et al. (2022) determined an 
upper limit for the isotropic prompt emission energy of 
AT2021lfa:; Eyiso < 1.2 x 10°! erg. The isotropic kinetic 
energy of AT20211fa, determined from the light curve fitting, is 
Exiso = 5.943 x 1074 erg. By applying the radiation efficiency 
Eso 


Ey iso + Ek iso 
radiation efficiency of AT2021lfa is 0.02%. Typically, 
observed GRBs display radiation efficiencies ranging from 
1% to 90%. However, the derived radiation efficiency for 
A'T20211fa falls significantly below this range, indicating a 
notably low efficiency. The exceptionally low radiation 
efficiency of AT20211fa may have two possible explanations. 
First, it may be an off-axis GRB whose prompt emission is not 
directed toward the Earth, so the observed radiation efficiency 
is very low. Alternatively, AT20211fa may be a “dirty fireball,” 
with its notably low radiation efficiency attributed to its low 
initial Lorentz factor. 


formula 7 = , we deduce that the upper limit for the 
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In order to avoid the so-called “compactness problem” and 
efficiently generate prompt gamma-ray emissions, the standard 
fireball model requires an initial Lorentz factor exceeding 100. 
Nevertheless, the circumburst environment around GRBs, no 
matter from the gravitational collapse of massive stars or the 
merger of compact binary stars, may carry an abundance of 
baryons, resulting in a challenge known as baryon contamina- 
tion. Under such circumstances, the abundance of baryons 
within the fireball will result in lower Lorentz factor known as a 
“dirty fireball.” While the prompt emission from “dirty 
fireballs” is faint and undetectable by GRB satellites, their 
multi-band afterglow remains observable, allowing us to 
directly detect their afterglow without relying on the prompt 
emission trigger. By utilizing relationships involving the 
Lorentz factor, peak timescale, interstellar medium, and 
isotropic equivalent kinetic energy of the jet (Nakar 2007), 
we can estimate the approximate value of the initial Lorentz 
factor as follows 

: ; 
5 8 
Hd (x + 2) : (2) 


Ib A zl 
tpeak = To 


Nism 


where Ex ;,o,5o represents the isotropic kinetic energy in units of 
10° erg, and z denotes the redshift. We have determined the 
burst time of AT202 I Ifa as Ty = 2021 May 3, 22:09:50 and the 
peak time as tpeak = 2021 May 4, 07:10:01. By substituting the 
fitted model parameters Ex iso,50 and n into Equation (2), we 
can estimate AT20211fa's initial Lorentz factor as about 18. 
This is consistent with the I'; — 20+ 10 derived by Lipunov 
et al. (2022), indicating that AT2021lfa may be a “dirty 
fireball,” and its low prompt emission efficiency may be 
attributed to its relatively small initial Lorentz factor. 
Additionally, the baryon loading of the fireball can be 
calculated using Ex is, = MI ge. Based on the bulk Lorentz 
factor distribution of GRB samples, Ghirlanda et al. (2018) 
found that the typical value of baryon loading is distributed 
around 10 5M... However, the baryon loading of AT20211fa is 
approximately 1.2 x 10 ^M. which is higher than the 
10 5M. obtained by Ghirlanda et al. (2018). This further 
supports the possibility that AT20211fa is a dirty fireball. 
Observations indicate that the GRB jets exhibit high 
collimation instead of isotropy. The prompt emissions of 
GRBs are more readily detectable only when the viewing angle 
(0,) is smaller than the opening angle of the jet core (0,) due to 
relativistic beaming effects. Consequently, a substantial portion 
of the prompt gamma-ray radiation from GRBs cannot be 
detected due to relatively large viewing angles. However, 
during the afterglow phase, as the Lorentz factor decreases, the 
radiation’s beaming angle widens, which enables the observa- 
tion of afterglow radiation to become possible even in the case 
of off-axis viewing. The study of off-axis afterglows assists in 
exploring the lateral structure of the jet (Ryan et al. 2015). 
Notably, GRB 170817A serves as an example of an off-axis 
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Figure 6. Correlations among timescales of GRB afterglow onset "bump." 


GRB, and multi-band afterglow light curve modeling suggests 
that its jet structure is more likely to be structured than uniform 
(Abbott et al. 2017; Alexander et al. 2018; Ghirlanda et al. 
2019; Lamb et al. 2019). In contrast, fiting the multi-band 
afterglow data of AT20211fa has revealed a viewing angle (0,) 
of 0.526 radians and a half-width of the jet core (6,) of 
0.662 radians. Therefore, it can be deduced that AT20211fa is 
observed on-axis. 

In the early phase of the GRB afterglow, interpreting several 
observed phenomena poses a challenge for the conventional 
external shock model because of the intricate interactions 
between the external shockwave and the circumburst medium. 
While conducting X-ray early afterglow fitting, Yost et al. 
(2003), Ioka et al. (2006), Fan & Piran (2006) and Wei & Fan 
(2007) postulated that the microphysical parameters e, and eg 
of the shockwave exhibited temporal variation rather than 
constancy. Subsequently, they found a significant improvement 
in fitting results by taking into account the time-dependent 
nature of these microphysical parameters. During the analysis 
of multi-band data for AT20211fa, it was observed that the late- 
time radio observations exhibit higher flux compared to the 
fitted light curve when keeping the microphysical parameters 
€», €g, and £y constant. To achieve a better fit to the 
observational data of AT20211fa, we have taken into account 
the temporal variation of these three parameters, resulting in a 
significant improvement in fitting the late-time radio observa- 
tions. The fitting results indicate a power-law decrease over 
time in the parameters €,, eg, and £y, which may be attributed to 
the gradual weakening of the shock's strength. 

Since the launch of the Swift satellite, an increasing number 
of early-phase afterglows have been observed. A smooth bump 
has been detected within the early light curves of some 
afterglows (Liang et al. 2010; Han et al. 2022). Liang et al. 
(2010) conducted a statistical analysis of afterglow samples 
exhibiting this bump feature and found that these bump 
structures exhibit similarities, indicating a strong correlation 
among the timescales of these bumps. The relationships among 


these timescales can be expressed as: 
log t; = (0.48 + 0.13) + (1.06 + 0.06)log t,, (3) 
logw = (0.05 + 0.27) + (1.16 + 0.10)log tp, (4) 


where w represents the full width at half maximum (FWHM) of 
the bump, t, denotes the time interval between the burst and the 
peak time, ¢, is the rising timescale measured at the FWHM, 
and fq is the declining timescale measured at the FWHM. Here 
we measure the characteristic timescales of the AT2021Ifa 
bump to be: w= 42,8328, t,=32,413s, t-=13,701s, and 
t4 = 29,131 s. We can compare the onset bump timescale of 
AT20211fa with that of the afterglow samples from Liang et al. 
(2010), as depicted in Figure 6. The timescale of AT20211fa 
adheres to the empirical timescale relationships obtained by 
Liang et al. (2010), but AT20211fa's bump displays a broader 
peak profile and a later peak time. This result is not surprising, 
as AT20211fa is highly likely to be a “dirty fireball,” leading to 
prolonged deceleration timescales and delayed peak times. 

Moreover, we conduct a comparison between the luminosity 
of AT2021lfa and the Swift GRB afterglows (Kann et al. 
2011), as depicted in Figure 7. To facilitate this comparison, we 
transformed the Swift GRB afterglow in the observer frame to a 
common redshift of z= 1. The observed luminosity of 
A'T20211fa falls within the range of observations for the Swift 
afterglow sample. However, at approximately 1 day after the 
burst, AT20211fa exhibits higher luminosity compared to the 
majority of Swift afterglow samples. This observation suggests 
that when searching for optical transient phenomena, ZTF 
prioritizes the detection of brighter sources. Consequently, 
studying faint optical afterglows requires the utilization of 
telescopes with higher sensitivity. 

A comparison of the physical parameters of AT2020blt, 
A'T2021any, and AT20211fa is presented in Table 2. The fitting 
results indicate that AT20211fa exhibits relatively large values 
for both the half-width of the jet core and the viewing angle. 
However, the observed GRBs generally display smaller half- 
widths of the jet core, along with smaller viewing angles. This 
trend could be attributed to observational selection effects. 
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Figure 7. Compare the brightness of Rc band between AT20211fa and Swift 
GRB samples. 


Table 2 

Comparison of AT2020blt, AT2021any, and AT20211fa Physical Parameters 
Transient AT2020blt AT2021any AT2021any AT20211fa 
0, [rad] 0.06007 0.557027 0.020.003 0.526 0181 
log; Eo 53.61*033 52.580083 52.9903. — 54774 40433 
6,[rad] 0.14004 0.96*031 0.11*001 0.662: 0297 
0, [rad] 0.05007 e 0.807045 z 

b 5.14529 " » * 
logiy to 1.902139 0.0651 -017 019 L039'0$9 
p 2,781914 2.30+905 2.35:005 3.092*0024 
logo € —-1407031 -1 -0.8:00j  —1.17749338 
logi €B — 1.641033 —2.23+013 —323 01$ — —4465'018 
ên 0.67703! 1 1 0.6957 9.528 
References Sarin et al. Gupta et al. Xu et al. this work 

(2022) (2022) (2023) 


Theoretically, scenarios should exist where both the half-width 
of the jet core and the viewing angle are relatively large (Gao & 
Dai 2010). Nevertheless, due to the nearly constant total energy 
released by the GRB central engine, an increase in the jet angle 
results in a decrease in the energy emitted per unit solid angle, 
which makes it more challenging to detect GRBs with larger jet 
angles. Our results indicate that AT20211fa is one of the optical 
transient sources with relatively large jet core half-widths and 
viewing angles. Furthermore, Gupta et al. (2022) adopted the 
final upper limit detection time as the burst time and found 
unusually large values for both the jet core half-width and 
viewing angle of AT2021any. In a recent study, Xu et al. 
(2023) employed the burst time as a free parameter to derive 
the jet core half-width and viewing angle values for 
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AT2021any, as presented in Table 2. The derived jet core 
half-width and viewing angle values for AT2021any closely 
resemble those obtained for AT2020blt, which indicates that 
the value of burst time is very important for describing the 
afterglow behavior accurately. 


5. Conclusions 


AT20211fa was detected by ZTF, followed by subsequent 
observations using Swift/XRT and ground-based telescopes. 
Notably, three hours before the initial ZTF observation, 
MASTER detected a gradual increase in the brightness of 
A'T20211fa. Due to its rapid flux decline in multiple wavelength 
bands, its reddish color (g — r0), and the absence of 
corresponding gamma-ray emission, AT2021lfa has been 
classified as a candidate for an “orphan afterglow.” 

In our study, we applied the empirical formula introduced by 
Kocevski & Liang (2001) to fit the onset bump of AT2021 Ifa, 
resulting in a burst time of Tọ — 2021 May 3 22:09:50 and a 
peak time of fg = 2021 May 4 07:10:01. Notably, when 
comparing our results to those obtained by Lipunov et al. (2022), 
our derived burst time is earlier, and the peak time is somewhat 
later. Subsequently, we employed afterglowpy to fit the 
afterglow observations of AT20211fa in the radio, optical, and 
X-ray bands after the peak time. During the fitting process, we 
considered power-law variations in microphysical parameters e,, 
€p, and £y over time, resulting in a significant improvement in 
fitting the late radio band data. We determined the half-width of 


the jet core 0. for AT20211fa to be 0.662*0297 radians and the 


viewing angle 0, to be 0,526104; radians, which does not 


support the possibility of off-axis observation for AT20211fa and 
suggests that AT2021lfa may be an on-axis source with 
relatively large jet and viewing angles. By using the relationship 
between the peak time and the initial Lorentz factor, we 
estimated AT2021Ifa’s initial Lorentz factor to be about 18, 
strongly suggesting that AT2021l1fa is likely a “dirty fireball.” 
Because of its relatively low initial Lorentz factor, the prompt 
emission is comparatively weak, making it undetectable by 
gamma-ray satellites. 

When comparing the initial bumps of AT20211fa with the 
afterglow sample, we found that they exhibit similar structures, 
although AT20211fa's bump is broader and reaches its peak at a 
later time. Additionally, we noted that AT20211fa's brightness 
falls within the range of brightness observed in the Swift 
afterglow sample. However, at about lday after the burst, 
AT20211fa's brightness exceeds that of most Swift afterglow 
samples, suggesting a bias toward observing brighter optical 
transient sources by the ZTF. In order to detect GRBs with 
weaker prompt emissions and dimmer afterglows, more 
sensitive high-energy satellites and optical telescopes, such as 
Space Variable Objects Monitor (SVOM) and the Einstein 
Probe, are needed to make significant contributions in this 
research field. 
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